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Abstract 

For injection of beams into circular machines with different horizontal and vertical emittance 
acceptance, the injection efficiency can be increased if these beams are flat, i.e. if they feature 
unequal transverse emittances. Generation of flat electron beams is well-known and has been 
demonstrated already in beam experiments. It was proposed also for ion beams that were generated 
in an Electron-Cyclotron-Resonance (ECR) source. We introduce an extension of the method to 
beams that underwent charge state stripping without requiring their generation inside an ECR 
source. Results from multi-particle simulations are presented to demonstrate the validity of the 
method. 

PACS numbers: 41.75.Ak, 41.85. Ct, 41.85.Ja 
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Beams provided by linacs generally have equal or quite similar horizontal and vertical 
emittances (e x and e y ). If beams from linacs are to be injected into a subsequent circular 
machine, the injection process might impose different acceptance limits (A x < A y ) on the 
two transverse emittances. One common injection scheme is the multi-turn injection (MTI) 
using a time dependent orbit bump in one plane for beam stacking in this plane. As a 
consequence the injected beam emittance in this plane, the horizontal for instance, must be 
significantly lower with respect to the other one. To meet the acceptance criteria in both 
planes, linac designers and operators try to keep both transverse emittances below the more 
stringent acceptance limit: 

e x -e y = e 2 x < A 2 X . (1) 

The requirement to the product of the emittances is relaxed if the linac beam will feature 
different transverse emittances as well, i.e. emittances being adopted to the acceptances of 
the machine into which injection has to be performed: 

A x < e x • t y < A x ■ A y . (2) 

To this end a single pass beam line is needed that provides transformation of a round 
uncorrelated beam (e x = e y ) to a flat beam (e x < e y ). Although for electrons round to flat 
transfer was proposed |l| and has been demonstrated experimentally j^J, to our knowledge 
it has not been proposed for initially uncorrelated round ion beams. Such an ion beam line 
must be different from an emittance swapping beam line that exchanges the two different 
transverse emittances by using three skew quadrupoles for instance 



I. REQUIRED MATHEMATICAL TOOLS 

A round to flat transformation implies beam line elements that cause coupling of the hor- 
izontal and vertical plane. It needs a four dimensional description. This section introduces 
the quantities referred to in the paper. The next section deals with a conceptual layout of a 
round to flat transformation for ions extracted from an ECR (Electron-Cyclotron- Resonance) 
source. Finally, a concept is presented that provides round to flat transformation for any 
uncorrelated ion beam being charge state stripped during transportation. GSI aims at ex- 
perimental verification of the concepts. Multi-particle simulations were done for sections 
planned to be integrated into the existing Universal Linear Accelerator (UNILAC) at GSI. 
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Any round to flat operation should only use beam line elements being linear in the four 
dimensional transverse phase space. In that case the transportation M of single particle 
coordinates from a position 1 to a position 2 is 
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mn "7,13 ™14 

m 2 i m 2 2 m 2 -i m 2 A 
m 3 i m 32 m 33 m 34 
m 4 i m 42 m 43 m 44 

As coordinates we use (x, y) to indicate transverse displacements and (x', y') to indicate 
their derivatives w.r.t. the position s along the design orbit. In the following sections we 
refer to beam dynamics simulation results from a code that also uses these coordinates. 
Additionally, the experimental proof-of-principles will be based on measurements of these 
coordinates. For this reasons we do not use the conjugate coordinates that additionally 
include the contribution of the magnetic vector field A in the momenta. 
The linear transport preserves the four dimensional rms emittance defined through the 
beam's second moment matrix 
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and 
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e 4d 



detC. 

Second moments of the beam are transported using the matrix equation 

(', = M C\ M T . 
Coupling between horizontal and vertical plane results in 



(5) 
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with equality just for zero inter-plane coupling moments. Here e x and e y are the rms emit- 
tances defined through the beam moments in their sub phase spaces in an analogue way to 
Equ. [5j Assuming an initial round beam without any inter-plane correlations, the trans- 
formation into a flat beam cannot be accomplished by applying symplectic transformations 



only [4|. Symplectic transformations M preserve 644, meet the condition 

M T JM = J 
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(9) 



and may just increase both transverse emittances (e x , e y ) while keeping them equal to each 
other. 

Rigorous analysis of coupled four dimensional transverse beam dynamics can be found in 
Q]. This section is concluded in providing eidetic justification for the beam line elements, 
which a round to flat transformation should include. In order to achieve emittance reduction 
in one plane, the horizontal for instance, a non-symplectic transformation M ns must be 
integrated into the round to flat transformation section. This transformation anyway should 
be linear to preserve e^. The section must also include an element M t that decreases the 
horizontal emittance and increases the vertical one. This element might act as 



5x = —a ■ x' 

5y' = b-y', (0,6) > 



(10) 



where a thin element not changing the particle position is assumed. The element imposes 
an effective linear friction in the horizontal plane and a linear heating in the vertical plane. 
If at the entrance to element M t the distribution has large inter-plane coupling moments 
< x'y > and < xy' > imposed by preceding elements, the approximation 



turns Eq. [TU] to 
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being the transformation of thin skew quadrupole with focusing strength 



< x' 2 > 
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(14) 



< xy' > < x' 2 > ' 
implying that < x'y > and < xy' > need to differ in sign. 

The required condition of different signs of the moments <x'y> and <xy'> at the skew 
quadrupole entrance can be achieved by including a transformation 
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through a solenoid fringe field into the beam line preceding the skew quadrupole with 
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B as the solenoid on- axis magnetic field strength, and (Bp) as the beam rigidity. The 
K values at the entrance and exit fringe are equal in magnitude but they differ in sign. 
Solenoid fringe fields can provide the non-symplectic transformation M ns that preserves 
e^d- Summarizing, a round to flat transformation of an initial beam without inter-plane 
correlations can be achieved if it includes solenoidal fringe fields causing a non-symplectic 
transformation M ns and skew quadrupoles for the emittance transfer. This does not mean 
that the transformation can be achieved using just two elements: a fringe field and a single 
thin skew quadrupole. Both elements are essential ingredients but they are not sufficient. 
Apart from this it is physically not possible to create a stand alone fringe field since static 
magnetic field lines are closed, i.e. an entrance fringe field intrinsically causes an exit fringe 
field. 



II. BEAMS EXTRACTED FROM AN ECR SOURCE 



Round to flat transformation of ion beams extractec 
same concept as this transformation of electron beams [l 



from an ECR source applies t 
,2]. For ions it was proposed in 



lie 



assuming a beam second moment matrix C with non-zero elements just along its two diag- 
onals. But generally all coupling moments are different from zero and their magnitudes are 
similar. For electron and for ion beam extraction the required non-symplectic transforma- 
tion M ns is provided through placing the beam generation inside a longitudinal field region. 
For electrons a solenoid is installed around the cathode. In case of ions the longitudinal 
field for plasma confinement is an intrinsic part of the ECR source. Prior to extraction the 
beam inside the plasma chamber is round (e x = e y ) and is not correlated. After extraction 
along the solenoidal fringe field the transverse phase space distribution from the source has 



inter-plane correlations. Fig. [T] shows such a distribution 



or 16 3+ extracted from the ECR 



source at GSI (Fig. [2]) with an energy of 2.5keV/u flQ. 
The distribution is obtained from multi-particle simulations and beam measurements per- 
formed at the HLI are in good agreement with these simulations 
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12|. In Fig. m (and 



in the subsequent ones of same type) correlations between two coordinates u and v are 
indicated by the dimensionless a u „-parameter defined through the moments as 

- <uv> 

V< u >< v 1 > — < uv > l 

The 120° symmetry in the inter-plane projections are due to the non-linear inter-plane 
coupling hexapolar field inside the ECR source. The observed inter-plane correlations a uv are 
from the solenoidal fringe field and the signs of the beam moments < x'y > and < xy' > differ. 
Although the absolute values of the inter-plane a-parameters seem small, the distribution is 
considerable correlated as seen by comparing the product of the two transverse emittances 
with the four dimensional rms emittance (see Eqs. [5] and [7]), i.e. e x ■ e y = 4.2 e^. Due to 
the intrinsic inter-plane correlations the product e x ■ e y exceeds the value required from four 
dimensional emittance preservation by a factor of about four. That implies in turn that 
decreasing e x ■ e y by this factor can be achieved by just removing the inter-plane correlations 
from the beam using linear beam line elements. 

Figure [3] displays a beam line that removes the inter-plane correlations from the beam. It 
comprises three solenoids, two quadrupoles, and two skew quadrupoles. Its total length 
is about 4 m, the focusing strength are less than 0.5 m _1 (normal and skew quadrupoles) 
and 4.5 m _1 (solenoids). After extraction from the source the beam is divergent in both 
transverse planes. To provide simultaneous focusing in both planes a solenoid is used as a 
first element of the beam line. The beam line needs to remove the four inter-plane correlation 
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FIG. 1. Two dimensional projections of the ECR transverse phase space distribution at the entrance 
to the emittance transfer section. Horizontal and vertical rms Twiss parameters are indicated as 
well as the four inter-plane correlation parameters. 
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IRON YOKE HEXAPOLE 



MICROWAVE 
14.5 GHz 




FIG. 2. (Color) ECR source at GSI with a permanent hexapole of 1 T at the pole tip. The 
longitudinal magnetic field is 1 T at maximum as well. It has an electron suppressing electrode at 
a potential of -2 kV. The extraction voltage is set to obtain 2.5 keV/u after extraction 11]. 



Solenoids Quadrupoles Skew Quadrupoles 

■ ■ D D ■ D I 

391 cm 

FIG. 3. (Color) Beam line from left to right for emittance transfer of a beam extracted from an 
ECR source (here the GSI HLI source). It comprises three solenoids, two quadrupoles, and two 
skew quadrupoles. Apertures are not to scale. 



moments, minimize the vertical emittance (for instance) at its exit, and provide for beam 
envelopes within a reasonable aperture. In order to meet these boundaries simultaneously, 
several beam line elements are needed and we found that the beam line in Fig. [3] can do so. 
Finding the required focusing strengths is accomplished by numerical tools based on linear 
transport of the second beam moments (Eq. [HD- As input the second moments extracted 
from the distribution shown in Fig. [T] has been used. Knowing the focusing strengths, multi- 
particle simulations using PARMTRA ]13) have been done. Finite magnets with hard edges 
were used. The energy spread of the beam is very low but has been considered anyway. 
Space charge is of no concern at the HLI ECR source. 
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FIG. 4. (Color) From top to bottom: horizontal beam envelope, vertical beam envelope, horizon- 
tal rms emittance, and vertical rms emittance along the emittance transfer section for the ECR 
distribution. 

Simulated beam envelopes together with transverse rms emittances along the beam line are 
plotted in Fig. HI The beam is round (equal emittances) until the exit of the third solenoid. 
At this location the distribution has maximum correlation with e x ■ e y = 12.5 e^. Despite 
this high value the beam is not subject to losses. The emittances are dominated by wide 
angular spreads. Within the final skew quadrupole the correlations are almost completely 
removed. The distribution at the exit of the section is shown in Fig. [5j Its vertical rms 
emittance is about 13 times smaller than the horizontal one. Compared to its initial value 
the vertical rms emittance was decreased by 86% while the horizontal one increased by 85%. 
The product of the two transverse emittances was reduced from 4.2 to 1.1 while 
was kept constant. The remaining factor of 1.1 as well as residual values of the inter-plane 
a-parameters shows that the correlations were not completely removed. We attribute this to 
non-linearities of the hexapolar coupling elements which cannot be fully removed by linear 
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FIG. 5. Two dimensional projections of the ECR transverse phase space distribution at the exit 
of the emittance transfer section. Horizontal and vertical rms Twiss parameters are indicated as 
well as the four inter-plane correlation parameters. 
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elements. Reducing correlations using hexapoles is proposed in 



III. BEAMS TO BE STRIPPED 

ECR sources generally do not provide high beam intensities for heavy and intermediate 
mass ions. Other source types are used in these cases (Multi-Cusp, Metal Vapor Vacuum 
Arc), which provide low charge states but high particle currents. From first principles also 
such sources can be designed applying ion extraction along a solenoidal fringe field to allow 
for later round to flat transformation. However, at extraction energy high intensity beam 
dynamics is subject to space charge forces. Latest at the entrance to an Radio Frequency 
Quadrupole (RFQ) eventual space charge neutralization does not mitigate any longer the ef- 
fect of beam self fields. Acceleration and transport with space charge might strongly reduce 
or remove the required inter-plane correlation that were imposed previously at extraction 
along a solenoidal fringe field. For this reason the beam should be exposed to the fringe 
field after some acceleration, when space charge forces are sufficiently reduced. On the other 
hand ion beams from low charge state, high particle current sources are stripped after gain- 
ing some energy in order to keep the acceleration process reasonable efficient. 
The obstacle to meet is to provide for a non-symplectic transformation M ns after acceler- 
ation. Unlike in the ion source, the accelerated beam must enter and leave the solenoidal 
field and a complete solenoid is a symplectic transformation which cannot be used to prepare 
round to flat transformation. 

The transformation through the solenoid is non-symplectic if the beam rigidity is abruptly 
changed in between the entrance and exit fringe fields, i.e. if the beam properties are re- 
set inside the solenoid. In this case the focusing strengths K from Eq. [15] are different at 
entrance and exit and the beam transport through the solenoid is a linear non-symplectic 
transformation M ns preserving e^. 

We propose to merge the required charge stripping of ions from low charge state, high par- 
ticle current sources and the provision of a transformation M ns for transverse emittance 
transfer preparation by placing the charge stripping foil in between two pancake coils as 
drawn in Fig. El Stripping abruptly changes the beam rigidity and has to be performed any- 
way along the ion linac. Based on the existing charge stripping section of the GSI UNILAC 
a beam line was conceptually designed that can provide horizontal to vertical emittance 
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Beam Axisz [arb. units] 



FIG. 6. (Color) Conceptual layout of an ion stripper set-up providing a non-symplectic transforma- 
tion M ns for later transverse emittance transfer. The stripping foil is placed between two pancake 
coils, i.e. inside a region with a longitudinal magnetic field of few Tesla in strength. 
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FIG. 7. (Color) The Universal Linear Accelerator (UNILAC) at GSI. 

transfer. 

The UNILAC (Fig. [7]) comprises three source terminals, two Drift Tube Linac (DTL) sec- 
tions separated by a gaseous stripper section, and a transfer channel to the synchrotron 
SIS18. Along the DTL sections the beam is accelerated to 1.4 MeV/u and 11.4 MeV/u, re- 
spectively. The transport channel to the synchrotron includes a foil stripper section shown 
in the photograph of Fig. [HJ The foil stripper section 14J comprises a quadrupole singulet, 
a stripping foil, a vertical four bend chicane, and a quadrupole doublet. The chicane pro- 
vides dispersion in its center for charge state separation. To demonstrate emittance transfer 
the beam line can be modified as drawn schematically in Fig. [9j GSI aims at experimental 
verification of the emittance transfer concept. The modifications of the existing foil stripper 
beam line will be conceived such that they do not intercept the beam in case no emittance 
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FIG. 8. (Color) Fish eye photograph of the foil charge state stripper at the GSI UNILAC including 
the four bend analyzing chicane. The beam enters from the left. 

Quadmpoles Chicane (4*35°) Skew Quadmpoles 

■ □ □ □ CLDCT 

Long. B-Field 
Strip. Foil (H 3 * -> 3p) 

1079 cm 

FIG. 9. (Color) Beam line from left to right for emittance transfer of a beam along the UNILAC. 
It comprises two quadrupole singulets, a stripping foil in between two pancake coils, a four bend 
chicane, a doublet, and a skew quadrupole triplet. Apertures are not to scale. 

transfer is required. The effective length of the longitudinal magnetic field region is 40 cm 
and the foil is placed 10 cm from the exit of this region. First beam tests are to be performed 
using a beam of Hjj~ stripped to a proton beam in a 20 fig /cm 2 carbon foil placed between 
to pancake coils. Protons were chosen in order to simplify the first proof of principles: 
low focusing strengths allow for usage of on-site equipment, low beam current avoids space 
charge effects, and a single charge state spectrum after stripping is provided. Multi-particle 
simulations for this case have been done and are presented in the following. 
The change of charge state, i.e. rigidity, inside a longitudinal magnetic field region is mod- 
eled by matrix transformation of single particle coordinates according to Eq. [3j The matrix 
includes to parts Msoh and Mg i o ; the first part comprising a solenoid entrance fringe field 
and the longitudinal magnetic field until the stripping foil: 
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with 

C := cos(KL), S := sin(KL), a := 2KL, (19) 

K from Eq. [TJ)J and L as distance from the effective longitudinal field edge to the foil. The 
beam rigidity is calculated from the unstripped charge state, i.e. using H3 at 11.4 MeV/u. 
The foil itself is modeled by increasing the spread of the angular distribution through scat- 
tering, i.e. 

< x' 2 > + < y' 2 > — ► < x' 2 > + < y' 2 > + Atp 2 . (20) 

The amount of scattering Aip was calculated using the ATIMA code [l^| as 0.296 mrad. 
Due to scattering the four dimensional emittance eu increases during the stripping process. 
Transport from the foil to the exit of the longitudinal magnetic field region is modeled by a 
matrix corresponding to Msok, i-e. 
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(21) 



using the beam rigidity corresponding to the stripped beam, i.e. protons at 11.4 MeV/u. 
Here L is the distance from the foil to the exit of the longitudinal magnetic field region. 
According to 15( the energy loss and straggling in the foil can be neglected. It must be 
emphasized that the increase of 644 is purely from scattering in the foil; it is not caused 
by the change of beam rigidity inside the longitudinal field. Scattering in the foil occurs 
anyway along beam lines that include foil strippers. The simulations additionally included 
beam momentum spread, hexapolar fringe fields of the chicane dipoles, and the existing 
aperture limitations. 

Simulations started at the entrance to the quadrupole singulet in front of the chicane. Fig. [TOl 
shows the initial transverse phase space distribution. Initial six dimensional Twiss param- 
eters were concluded from various beam experiments with Hjj" beams. Beam spot images 
from fluorescence screens close this position did not reveal any x-y correlation. Therefore 
initial inter-plane correlations are assumed as zero. The distribution was tracked through 
the emittance transfer section and corresponding envelopes and transverse emittances are 
plotted in Fig. [TTJ Full beam transmission is achieved. The observed increase in both 
transverse emittances at the entrance to the longitudinal magnetic field region is driven by 
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FIG. 10. Two dimensional projections of the transverse phase space distribution at the entrance 
to the proposed emittance transfer section along the UNILAC. Horizontal and vertical rms Twiss 
parameters are indicated as well as the four inter-plane correlation parameters. 
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FIG. 11. (Color) From top to bottom: horizontal beam envelope, vertical beam envelope, horizontal 
rms emittance, and vertical rms emittance along the proposed emittance transfer section along the 
UNILAC. 



inter-plane correlations caused by the fringe field. There is no increase of the four dimen- 
sional emittance until the stripping foil. However, the increase of at the foil from 
scattering is just 4%. Vertical dispersion from the first dipole causes additional increase of 
the vertical and the four dimensional emittance, both being compensated later by the last 
dipole. 

The distribution at the exit of the longitudinal magnetic field region is shown in Fig. [121 It 
has inter-plane correlations and the signs of < x'y > and < xy' > differ as expected from the 
action of the solenoidal fringe fields. The product of the transverse emittances is a factor of 
1.8 larger than the four dimensional emittance e^. During the transport through the chicane 
the horizontal emittance is preserved. The intrinsic inter-plane correlation of the beam, i.e. 
the ratio e^y/e^ (corrected w.r.t. vertical dispersion) remains constant. Along the elements 
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FIG. 12. Two dimensional projections of the transverse phase space distribution at the exit of 
the longitudinal magnetic field section of the proposed emittance transfer section along the UNI- 
LAC. Horizontal and vertical rms Twiss parameters are indicated as well as the four inter-plane 
correlation parameters. 
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after the chicane the correlations are removed, the horizontal emittance is minimized, and 
the beam is re-matched for further transport to the synchrotron. To accomplish these tasks 
within the existing apertures, the gradients of the quadrupole doublet, the skew triplet, and 
the final singulet are determined numerically. The same algorithm as mentioned in Sec. [TT] 
is applied. As input for the algorithm the second beam moments extracted from the distri- 
bution at the doublet entrance are used (Fig. [13]). Here we presently rely on the simulations 
using the PARMTRA code being part of the PARMILA code family The later has 

been benchmarked successfully with beam experiments involving more complex beam lines 



even with strong space charge 







181 ] . The beam moments at the doublet entrance also 



can be measured if a pepper-pot device will be installed right behind the last singulet. Its 
installation is planned for the experimental verification of the proposed emittance transfer 
concept. 

After applying the required gradients for emittance transfer and horizontal emittance mini- 
mization the beam is transported until the exit of the last singulet. The resulting distribution 
is plotted in Fig. [TU Inter-plane correlations are almost completely removed. Residual val- 
ues are due to hexapolar dipole fringe fields as well as to chromaticity in quadrupoles and 
skew quadrupoles from the finite beam momentum spread (~ 0.1% rms). The product of 
e x ■ e y exceeds by less than one percent. Starting from an initial transverse emittance 
ratio of 1.0, the final ratio is 4.8. The horizontal emittance is reduced by 53% and the 
vertical one increased by 125%. The four dimensional emittance increased by 5% mainly 
due to angular scattering in the stripper foil (4%). The remaining part is from beam line 
non-linearities as mentioned before. 

Unlike for the round to flat transformation of beams from an ECR source (Sec. [TTJ) , the 
longitudinal field strength at the stripper foil is a free parameter. The amount of emittance 
transfer scales with the field strength. However, practical limitations for this amount rise 
from finite apertures. 



IV. CONCLUSION AND OUTLOOK 

It was shown that flat to round transformation is not restricted to electron beams but 
can be accomplished for beams of ions as well. The required solenoid fringe fields are 
provided intrinsically for beams being extracted from an ECR source. Beams from other 
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FIG. 13. Two dimensional projections of the transverse phase space distribution at the entrance 
to the first quadrupole after the chicane of the proposed emittance transfer section along the 
UNILAC. Horizontal and vertical rms Twiss parameters are indicated as well as the four inter- 
plane correlation parameters. 



19 



-i 1 1 r- 



5 - 



T3 



e = 0.56 mm mrad 
-5 hp= 14.6 m 
a = -1.66 



_i i i i_ 



-20 20 

x [mm] 



20 



E 
E 



-20 - 



- 1 1 1 1 r~ 



a = -0.05 

. I I I I I L_ 



s = 1079 cm 



-20 
x [mm] 



20 



T3 
03 



E 




-20 
x [mm] 



T3 



E 




e 4d = 1 .51 mm mrad 



20 
2 y [mm] 



T3 




x' [mrad] 



T3 



- 1 — i — i — i — i - 




e = 2.71 mm mrad 
B= 8.47 m 
a = -1.01 



-20 20 
y [mm] 



FIG. 14. Two dimensional projections of the transverse phase space distribution at the exit of 
the proposed emittance transfer section along the UNILAC. Horizontal and vertical rms Twiss 
parameters are indicated as well as the four inter-plane correlation parameters. 
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sources can be emittance-shaped by applying charge state stripping inside a longitudinal 
magnetic field. Examples for the layout of emittance transfer beam lines were presented for 
both cases. Multi-particle simulations demonstrated that the transfer is feasible. An initial 
round beam could be transformed into a beam with transverse emittance ratio of almost 
five. The amount of transfer for stripped ions can be controlled through the longitudinal 
magnetic field strength. 

GSI aims at experimental verification of emittance transfer by assembling the beam lines 
presented in this paper. Emittance transfer of beams from ECR sources might by of special 
interest for facilities operating an ECR source in connection with multi-turn injection into 
a synchrotron, as cancer treatment facilities for instance 19] . High intensity beams of i 
mediate or heavy mass ions for synchrotrons are to be provided for the FAIR project 
Stripping of those ions results in a charge state spectrum and requires charge state sepa- 
ration. Additionally, stripping of heavy ions uses much thicker foils (~ 600 ng/cm 2 ) with 
respect to light ions, thus the four dimensional emittance growth from scattering is sig- 
nificantly higher (> 70%). This growth reduces the budget for the amount of emittance 
transfer. However, emittance transfer lines can be designed also for such scenarios although 
the layout will be more complex than the one for protons presented here. Tentative layouts 
were made on the emittance transfer for beams of 238 U 27+ stripped to 238 U 73+ based on the 
existing foil stripping section at the GSI UNILAC. Horizontal emittance reduction of up 
to 40% was achieved including space charge and charge state separation. Further optimiza- 
tion of the layout for an emittance transfer section for intense beams of heavy ions is needed. 
This optimization is beyond the scope of the present paper. 
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